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Abstract 

We studied properties of the strange axial mesons in the relativized quark 
model. We calculated the Ki decay constant in the quark model and showed 
how it can be used to extract the Kii^Pi) — Ki(^Pi) mixing angle {Ok) 
from the weak decay r — > KiVt-. The ratio BR{t Vt-Ki{\21{))) / BR{t — > 
i/T-i^i(1400)) is the most sensitive measurement and also the most reliable 
since the largest of the theoretical uncertainties factor out. However the cur- 
rent bounds extracted from the TPC/Two-Gamma collaboration measure- 
ments are rather weak: we typically obtain —30° < Ok ~ 50" at 68% C.L. We 
also calculated the strong OZI-allowed decays in the pseudoscalar emission 
model and the flux-tube breaking model and extracted a ^Pi Pi mixing 
angle of Ok — 45". Our analysis also indicates that the heavy quark limit 

does not give a good description of the strange mesons. 
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I. INTRODUCTION 



The strange axial mesons offer interesting possibilities for the study of QCD in the non- 
perturbative regime through the mixing of the ^Pi and ^Pi states. In the SU(3) limit these 
states do not mix, just as the ai and bi mesons do not mix. For a strange quark mass greater 
than the up and down quark masses SU(3) is broken so that the ^Pi and ^Pi states mix 
to give the physical Ki states. In the heavy quark limit where the strange quark becomes 
infinitely heavy, the light quark's spin couples with the orbital angular momentum resulting 
in the light quark having total angular momentum j = | in one state and j = ^ in the other 
state, each state having distinct properties By studying the strange axial mesons 

and comparing them to the heavy quark limit one might gain some insights about hadronic 
properties in the soft QCD regime. 

Recently, the TPC/Two-Gamma Collaboration has presented measurements for the de- 
cays T- z/^Kf (1270) and z/^iTf (1400) [|]. It is expected that the LEP, CLEO, and 
BES collaborations, with their large samples of r's, will be able to study these decays in 
further detail These decays provide another means of studying ^Pi Pi mixing of the 
strange axial mesons in addition to using their partial decay widths and masses. 

In this paper we study the properties of the strange axial mesons in the context of the 
relativized quark model We compare the experimental measurements to the predictions 
of the model to extract the ^Pi — ^ Pi mixing angle (Ok)- Comparing both the experimental 
measurements and model results to various limits helps in understanding the nature of QCD 
in the soft regime. 

We begin in Sec. II with a brief description of the relativized quark model and a de- 
scription of the ^Pi Pi mixing. By comparing the mass predictions of the quark model 
to the observed Ki masses we obtain our first estimate for 6k- In Sec. Ill we calculate the 
Ki decay constants using the mock-meson approach and use the results to obtain a second 
estimate of 9k- In Sec. IV we study the strong decay properties of these states using the 
pseudoscalar emission model and the flux-tube breaking model Q and use the results 
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as another way of measuring the ^Pi Pi mixing angle. When appropriate we examine 
the non-relativistic and heavy quark hmits to gain insights into the underlying dynamics. 
Various aspects of the phenomenology of the strange axial mesons have also been studied 
by Suzuki in a series of recent papers P,p^ using approaches complementary to ours. 



II. THE Ki MASSES AND ^Pi Pi MIXING 

In this section we give a very brief description of the relativized quark model 
The spin-orbit contributions in particular will be important in understanding the ^Pi Pi 
mixing. The model is not derived from first principles but rather is motivated by expected 
relativistic properties. Although progress is being made using more rigorous approaches, the 
relativized quark model describes the properties of hadrons reasonably well and presents an 
approach which can give insights into the underlying dynamics that can be obscured in the 
more rigorous approaches. 

The basic equation of the model is the rest frame Schrodinger-type equation. The ef- 
fective potential, Vqq{p,r), is described by a Lorentz- vector one-gluon-exchange interaction 
at short distances and a Lorentz-scalar linear confining interaction. Vqg{p, r) was found by 
equating the scattering amphtude of free quarks, using a scattering kernel with the desired 



Dirac structure, with the effects between bound quarks inside a hadron [|Tl|]. Due to relativis- 
tic effects the potential is momentum dependent in addition to being co-ordinate dependent. 
The details of the model can be found in Ref. [||. To first order in (f/c)^, Vqq{p,r) reduces 
to the standard non-relativistic result: 

Vqq{p, r) V{r) = ff-"/ + if-"* + Hlf + H^^' (1) 

where 

H^l^f = C + br + ^Fq ■ Fq (2) 
includes the spin-independent linear confinement and Coulomb-like interaction. 
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the Thomas precession term. In these formulae, {Fq ■ Fg) = —4/3 for a meson and ^^(r) is 
the running coupling constant of QCD. 

For the case of a quark and antiquark of unequal mass the ^Pi and ^Pi states can mix 
via the spin orbit interaction or some other mechanism. Consequently, the physical j = 1 
states are linear combinations of ^Pi and ^Pi which we describe by the following mixing: 



K+ 



(^Pi)+ COS 9 + (2Pi)+ sin e = K+ cos 9 + K+ sin 9 
KtiHqh) = - ('Pi)+ sin ^ + ('P,y cos 9 = -X+ sin 9 + cos 9 



(8) 



The Hamiltonian problem was solved using the following parameters: the slope of the 
linear confining potential is 0.18 GeV^, m„ = = 0.22 GeV and = 0.419 GeV. The 
resulting masses of the unmixed states are: 



M{Ka) = 1.37 GeV 
M{Kk) = 1.35 GeV. 



(9) 
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We expect these values to be reasonable estimates as the model's predictions for the closely 
related ai and 61 masses are consistent with the experimental measurements. In this model 
spin-orbit mixing results in Ok = —5" but the Ki masses remain the same within the 
given numerical precision. These mixed masses and the mixing angle are not consistent with 
the measured values. 

We can obtain a phenomenological estimate of 9k by considering the 2x2 matrix relating 
Ka and Ki, to the physical i^i's. We do not make any assumptions about the origin of the 
^Pi Pi mixing and treat the off-diagonal matrix element of the Ki mass matrix as a free 
parameter. Diagonalizing the Ka — Kb mass matrix gives the relation between 6k and the 
mass differences: 

Oft M(ir,) - MjKb) 

coszUk - M(KiiU02)) - M(iri(1273)) ^ ^ 

with corresponding Ki masses: 

sin^ 2eK 



Mi,^ = Mb cos" Ok + sin' Ok - (M„ - M5 



2cos2^;^ 



sin*^ 2,9 

M^gh = Mb sin^ Ok + M, cos^ Ok + (M, - Mb)- (11) 

2 cos 26 k 

Solving gives 9k — ±41°. Note that degenerate Ka and Kb masses will always result in a 
mixing angle of ±45° [|12| . Thus, the value we obtain for 6k is more a reflection of the near 
degeneracy of the model's prediction for Mk^ and Mk^ than anything else and one should 
not read too much into the value we extract here. 



III. WEAK COUPLINGS OF THE Ki'S 

We use the mock meson approach to calculate the hadronic matrix elements ||6|jl3|-[17| . 
The basic assumption of the mock meson approach is that physical hadronic amplitudes 
can be identified with the corresponding quark model amplitudes in the weak binding limit 
of the valence quark approximation. This correspondence is exact only in the limit of zero 
binding and in the hadron rest frame. Away from this limit the amplitudes are not in 
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general Lorentz invariant by terms of order Pi/m^. In this approach the mock meson, which 
we denote by M, is defined as a state of a free quark and antiquark with the wave function 
of the physical meson, M: 

\M{K)) = y2M~ J d^p <^M{p)Xss<Pgg<Pcoiour\q[{mg/ii)K + p, s] q[{mg/ii)K - p, s]) (12) 
where $m(p), Xss, 4'qq-i ^^'^ 4>coiour are momentum, spin, flavour, and colour wave functions 

— * 

respectively, fj, — ruq + m^, K is the mock meson momentum, is the mock meson 
mass, and ^2M^ is included to normalize the mock meson wavefunction. To calculate the 
hadronic amplitude, the physical matrix element is expressed in terms of Lorentz covariants 
with Lorentz scalar coefficients A. In the simple cases when the mock- meson matrix element 
has the same form as the physical meson amplitude we simply take A ~ A. 
In the case of interest, the axial meson decay constants are expressed as: 

- l,)q I M(i?, A)) = A) (13) 

where e^{K, A) is the Ki polarization vector and fx^ is the appropriate Ki decay constant. 
To calculate the left hand side of Eq. (13) we first calculate 

(0|g7'^(l-75)g \q[{mJi^)K + p,s] q[{mg/ fi)K - p,s]) (14) 

using free quark and antiquark wavefunctions and weight the result with the meson's mo- 
mentum space wavefunction. 

There are a number of ambiguities in the mock-meson approach and different prescrip- 
tions have appeared in the literature. For example, there are several different definitions for 
the mock-meson mass {M^) appearing in Eq. 12. To be consistent with the mock meson 
prescription, we should use the mock meson mass defined as (Eg) + {Eg). However, because 
it is introduced to give the correct relativistic normalization of the meson's wavefunction 
the physical mass is another, perhaps more appropriate, definition. The second ambiguity 
is the question of which component of the 4- vector in Eq. 13 we should use to obtain fx^- 
In principle, it should not matter as both the left and right sides of Eq. 13 are Lorentz 
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4-vectors. This is true in the weak binding hmit where binding effects are totally neglected, 
but in practice, this is not the case. We follow Ref. [0 and extract fxi using the spatial 
components of Eq. 13 in the limit K 0. Finally, evaluating Eq. 14 introduces factors 
of ~ iTLi/Ei. While some prescriptions take the expression derived from Eq. 13 only as a 
guideline and introduce powers of ~ {mi/EiY with e an arbitrary power, we chose to use the 
expression exactly as derived from take Eq. 13. The different prescriptions are described 
in greater detail in Ref. [|l3l which calculated the pseudoscalar decay constants (fx)- We 
will follow the approach taken there and use the variations in prescriptions as a measure 
of how seriously we should take our results. In our results we therefore use the "exact" 
expression for f^^ and we take M^^ to be equal to the physical mass (Mphys). Variations in 
the mock-meson normalization result in variations in fx^ of at most 20%. Results using the 
physical mass lie in the middle of the range so that we expect uncertainties introduced by 



taking M^g = Mphys to be no more than ~ 10%. As in Ref. [|T3| fx^ was most sensitive to 
the wavefunction used. Here we use the sets of wavefunctions that gave the best agreement 
with experiment for in Ref. |TB|. We choose two possibilities, one which underestimated 
fx and one which overestimated it. We would expect these choices to likewise bound the 
actual value of the /^i. 

The expressions we obtain for fx^ are given by: 



4v^ 



M 



ji^2 [^^<PkM \ P 



Ea + ma\ I Eg + mg 



2E, 



2Eg 



1/2 r 



+ 



Eg + mq Eq + rriq 



Ki 



(2^ 1 P 



X 



Eq + mq\ f Eq + niq 



2Er, 



2Eg 



1/2 



1 



Ea + ma Eg + Ulg 



(15) 



where (j)Ki{p) is the radial part of the momentum space wavefunction, Eq = J\p\'^ + m'^ and 



E, 



\p\'^ + m?. In the SU{3) limit only Ka couples to the weak current. 
With the definition of fxi given by Eq. (13) the partial width for r KiUr is given by: 



A. The Non-relativistic Limit 



It is useful to examine the Ki decay constants in the non-relativistic limit where their 
qualitative properties are more transparent. In this limit the axial- vector meson decay 
constants become: 
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(17) 



r=0 



where Rp{r) is the radial part of the coordinate space wavefunction. Combining the weak 
decay amphtudes with the mixed Ki eigenstates the decay constants for the mixed \K^) = 
— \us) states are given by:[] 
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where we have defined 



A 



6M 



3 dRp{r) 



^1 I V Svr dr 



(19) 



r=0/ 



In the SU(3) limit /xi(^-Pi) explicitly goes to zero and only the ^Pi state couples to the 
weak current. The Kb coupling therefore goes like the SU{3) breaking {nis — md). 



^Note that the sign change going from Eqn. 17 to Eqn. 18 comes from the phase in the 
flavour wavefunction. 
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B. Extracting 9k Using the Non-Relativistic Expressions 



We can obtain an estimate of the ^Pi Pi mixing angle by comparing the quark model 
predictions to experiment. As stated above, the values of the decay constants were quite 
sensitive to the choice of wavefunction. We calculated the fx^ for two sets of wavefunctions 
that gave the best agreement between a quark model calculation and experiment for the 
pseudoscalar decay constants [O. We expect that the actual values for the fx^ will lie 



between the values predicted using these wavefunctions. The values for the two meson 
masses and two sets of wavefunctions are given in Table I. 

There are four measurements that can be used to constrain 9k- The TPC/Two-Gamma 
Collaboration 0] has made the measurements: 

BR{t vKi{mQ)) = (0.411;^^) x 10"^ (20) 

BR{t uKiiUOO)) = (O.TQtil) x 10"^ (21) 

BR{t vKi) = (1.17^:^^) X 10"^ (22) 



and Alemany |jT9|] combines CLEO and ALEPH data to obtain: 

BR{t uKi) = (0.77 ± 0.12) x 10"^ (23) 

which is smaller than, but consistent with, the TPC/Two-Gamma result. CLEO claims that 
the r decays preferentially to the 7^^1(1270). 

Using the ratio BR{t — > vKi{127Qi)) / BR{t z/fCi(1400)) has the advantage of factoring 
out the uncertainties associated with the Ki wavefunction. The ratio is given by0 



R = 1.83 



2 

sin 6k — S cos 6k 



cos 6k + S sin 6k 

where 1.83 is a phase space factor and 6 is an SU{3) breaking factor given by 



(24) 



^The numbers from Table I give a slightly different value since the different Ki masses in our 
expression for /ki do not exactly factor out. 
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The ratio, R, is plotted in Fig. la as a function of 9k- Taking ruu = 0.33 GeV and = 
0.55 GeV and fitting Eq. (24) to the ratio of the TPC/Two-Gamma resuhs we obtain 
—30° ^ 9k ^ 50° at 68 % C.L. where the large uncertainty is directly attributed to the large 
errors in the branching ratios.0 

Although the relative errors for the individual branching ratios are smaller than those 
of the ratio, especially for the sum to the two Ki states, using the branching ratios intro- 
duces additional uncertainties due to the errors associated with the meson wavefunction. 
In addition, the branching ratios turn out to be less sensitive to 9^ than the ratio. This 
is seen very clearly in Figs, lb, Ic, and Id where we have plotted the branching ratios 
for r z/T-iri(1270), r UrKi{lAOO) and the sum of the two respectively. The values 
= (295.6 ± 3.1) x 10"^^ s and |Ks| = 0.2205 ± 0.0018 were used to obtain these curves 
2T| . The two curves in each figure represent the two wavefunctions we use and we have 



included the experimental value with its error. In Fig. Id both the TPC/Two-Gamma and 
the CLEO / ALEPH values are shown. It is apparent from these figures that it is not partic- 
ularly meaningful to extract a value for 9k from these results and any value would be very 
model dependent. Clearly better data is needed. The ratio of the rates into the individual 
final states will give the most model independent constraints on 9k- 



C. Extracting 9k Using the Relativized Expressions 

We next calculate the axial meson decay constants using the relativized formula of 
Eq. (15). One might question the importance of including relativistic corrections. However, 



^ The of the fit actually has 2 local minima corresponding to both a negative and positive 
solution. However, since the hump separating the two solutions is approximately equal to Ax^ < 1 
the entire range given for 6k is consistent at 68 % C.L. 
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we need only consider the importance of another relativistic correction: QCD hyperfine in- 
teractions which give rise to the p — vr, K* — K, . . ., B* — B sphttings Although it is 
difficult to gauge the importance of relativistic corrections to the fx^, if nothing else their 
inclusion acts as one more means of judging the reliability of the results. 

As in the previous section we give results for two wavefunction sets that give reasonable 
agreement for the J'k in a similar calculation. The various Jki are given in Table II. We 
expect that the actual values will lie between the two values given for each case. The 
predictions for the various branching fractions are shown in Fig. 2 as a function of 6k along 
with the experimental values. The most reliable constraint again comes from the ratio of 
branching fractions which gives —35° ^ 9k ^ 45° at 68 % C.L. One could also extract 
values using BR{r z/_ft'i(1400)) and BR{t — > z/_ft'i(1270)) but as in the non-relativistic 
case these values are quite sensitive to the magnitude of which depends on the poorly 
known Ki wavef unctions. 

We conclude that the decays r vKi offer a means of measuring the ^Pi Pi mixing 
angle but to do so will require more precise measurements than are currently available. 



IV. STRONG DECAYS OF THE Ki'S 

It is well known that the strong decays of the Ki mesons provides a means of extracting 
the -1 Pi mixing angle [||]. In particular the BR{Ki{1270) K*n)/ BR{Ki{UOO) 
K*7i) and BR{Ki{l270) Kp)/BR{Ki{im) -> Kp) ratios have been especially useful. 
We examine the decays to the final states Kp, Ku, and K*7[. Although other decays 
are observed they lie below threshold and proceed through the tails of the Breit-Wigner 
resonances making the calculations less reliable. In this section we examine the strong Ki 
decays using the pseudoscalar emission model 0, the ^Po model (also known as the quark- 
pair creation model) [Q, and the fiux-tube breaking model |p. We concentrate on the 
decays Ki —>■ K*tt and Ki —>■ pK. 

For the decays Ki — > VP, where V and P denote vector and pseudoscalar mesons 
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respectively, the OZI-rule- allowed decays can be described by two independent S and D- 
wave amplitudes which we label S and D. The decay amplitudes, using the conventions of 
Eq. (8) are given by 



^ [K*7t]s) = -SsiniOK - 



0) 



A{K'r [K*7:]o) = +Dcos{9k - 9o) 



. 0) 

A{K'r ^ [pK]s) = +SsiYi{eK + ^o) 
A{K'r ^ [pK]d) = +Dcos{9k + ^o) 
A{K^^^^ ^ [pK]s) = +Sco^{dK + ^^o) 

A{K['''' ^ [ujK]d) = -JId sin(^^;, + ^o) (26) 



where sin^o = y 1/3 and cos^o = y2/3 and the subscripts S and D refer to S- and D-wave 
decays. In the heavy quark limit the j = 1/2 state decays into K*Tr in an S-wave and the 
j = 3/2 state decays into K*tt in a D-wave. Since the decay -ft'i(1400) — > i^*7r is dominantly 
S'-wave while the decay fi'i(1270) K*tx has comparable S and D-wave contributions we 
conclude that experimental data favours the heavier = 1+ to be mainly j = 1/2 and the 
lighter one to be mainly j = 3/2. 

In the following sections we give results for these amplitudes, the resulting decay widths 
and the fitted values of 6k for the various decay models. 



A. Decays by the Pseudoscalar-Meson Emission Model 



In this approach meson decay proceeds through a single-quark transition via the emission 
of a pseudoscalar meson 0. We assume that the pair creation of u, d, and s quarks is 

12 



approximately SU{3) symmetric. We follow Ref. [^] and use the various approximations 
introduced there. The resulting amplitudes are given by 



where A = 1.67, S = 3.27, q is the momentum of each outgoing meson in the centre of mass 
(CM) frame, q = q/f3, f3 = 0.4 GeV and 



Numerical values for the amplitudes are given in Table III. 

The partial widths for Ki K*tt and Ki —>■ pK and the ratio of the D to S amplitudes 
for Ki K*TT are plotted in Fig. 3 as a function of Ok for the i\:i(1270) and is:i(1400). The 
experimental values are given with their errors. From the figures it is clear that the experi- 
mental values for Ki(1270) ^ K*7c, Ki{UOO) pK, and Ad/As[Ki{UOO)] correspond to 
minima in the quark model results with 6^ ~ 45°. We performed a fit to the data listed 
in Table IV and obtained 9k = 48° ± 5°. We also allowed the S, A, and /5 parameters to 
vary and obtained very similar results, the main difference being that the value at the 
minimum decreased significantly. The partial widths and A£,/As{Ki — > K*7t) ratios are 
given in Table IV for the fitted value of 6k- 



The flux-tube breaking model is a variation of the ^Pq niodel which more closely describes 
the actual decay processes. In the ^Pq model the elementary process is described by the 
creation of a qq pair with the quantum numbers of the vacuum, J^*^ = 0^+. The greatest 
advantage of this approach is that it requires only one overall normalization constant for the 
pair creation process. In the fiux-tube breaking model, the fiux-tube-like structure of the 
decaying meson and its implications for ^Pq amplitudes are taken into account by viewing a 



D = -AfFiq') 
S = SF{q^) 



(28) 



(27) 




(29) 



B. Decays by the Flux- Tube Breaking Model 
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meson decay as occurring via the breaking of the flux-tube with the simultaneous creation of 
a quark-antiquark pair. To incorporate this into the ^Pq model, the pair creation amplitude 
7 is allowed to vary in space so that the qq pair is produced within the confines of a flux- 
tube-like region surrounding the initial quark and antiquark. This model is described in 
detail in Ref. The ^Pq model corresponds to the limit in which 7 is constant. 

For the ^Pq model using simple harmonic oscillator wavefunctions the S and D amplitudes 
are given by: 



S 

D = q 



^11,(3% + 17121,(31; ' 

2 i ^nP's + ^23(^0 
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F{q')A 



3(31,(3^ 
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277^1/4^1/2 ^(3^) ^(3^(3^ 
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mi + ma 
m2 



(32) 

(33) 

(34) 
(35) 
(36) 



m2 + ma 

mi and m2 are the quark and antiquark masses from the original meson, ma is the mass of the 
created quark/antiquark, the I3i are the simple harmonic oscillator wavefunction parameters, 
and q is the momentum of each outgoing meson in the CM frame. For these results we take 
the Mi to be equal to the calculated masses of the mesons in a spin-independent potential 
P]. Numerical values for the relevant amplitudes are given in Table III. 

The decay amplitudes in the ^Pq model were computed symbolically using Mathematica 



25|. In the flux-tube breaking model two of the six integrals were done analytically; the 



remaining four were done numerically. The integrands were prepared symbolically using 
Mathematica and then integrated numerically using either adaptive Monte Carlo (VEGAS 
2^) or a combination of adaptive Gaussian quadrature routines. 
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We calculated the Ki strong decays using both the flux-tube breaking model and the 
^Pq model for several sets of wavefunctions. In all cases we fitted 7 to 28 of the best known 
meson decays by minimizing the defined by = Y.ii^i'^-VTy/STj where 5T, is the 
experimental error The details of these fits are given in Ref . |^ . We performed a second 
fit to the Ki decays where we allowed both 9k and 7 to vary. The value of 9^ obtained in 
the second approach did not change much from the first value — the main difference was 
that the in the second fit was reduced substantially. The values for 7 obtained in the 
second set of fits are consistent, within errors, with those obtained by the global fit of Ref. 



27| . In Fig. 4 we show the decay widths and ratios of D to S amplitudes as a function of 
9k for the ^Pq model. The results for the two variations of the fiux-tube breaking model 
are very similar and are therefore not shown. It is clear from these figures that 9k will be 
approximately equal to 45°. The fitted values of 9k for the various models, and the resulting 
widths, are given in Table IV. 



V. DISCUSSION 



One of the motivations for this analysis is to relate hadron properties to the underlying 
theory via effective interquark interactions |^ . We begin our discussion of the Ki mesons 
by rewriting the non-relativistic spin dependent potential in a more suitable form and in- 
terpreting it as an effective interaction pSf . We will later examine the Ki meson properties 
in the limit tuq 00. 

The spin-orbit Hamiltonian can be rewritten as: 

2 
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(37) 
(38) 



^For the calculations in the flux-tube breaking model, a 1% error due to the numerical integration 
was added in quadrature with the experimental error. 
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where S = Sg + Sg, = Sg — Sg. Taking q = Q the various terms in H^ °- can be rearranged 
as 

= S ■ L + H^^ ■ L 

= Hl^, Sg-L + i/Q . Sq ■ L (39) 

where the definitions of H^^ , H^^ , H^^ , and Hf^ follow from Eqs. (37) and (38). It is the 
term which gives rise to the spin-orbit mixing between the singlet and triplet states. 
With this Hamiltonian, we obtain the following mass formulae for the P-wave mesons: 

•M(3Pi)\ f Mo + \{H,ont) + ^{Hten) - {Hto.) V^iH^J \ ( ^ 

M(3Po) = Mo + \{H,ont) - {Hten) " 2{H^J (40) 

where the (Hi) are the expectation values of the spatial parts of the various terms, Mq is 
the center of mass of the multiplet, and we have adopted a phase convention corresponding 
to the order of coupling L x Sg x Sq. 

We can rewrite H^°- using the substitutions rh = ^{nig + mq) and A = {niQ — nig) to 
obtain the approximate expression 

A 



8^ l_ldH^ 

3 m^r^ 2m^ r dr 



S-L- _ 
m 



^la^n 

3 m?r^ 2fn^ r dr 



S. ■ L. (41) 



Written in this way one sees that there is a factor of two difference between the colour 
magnetic term and the Thomas precession term for the relative to H'^^. The observed 
spin-orbit splittings in hadrons indicate a delicate cancellation between the colour magnetic 
and Thomas precession spin-orbit terms. Given this cancellation, the factor of two could 
lead to a large effect or even a sign reversal in the spin-orbit mixing. 

In particular, the relativized quark model gives Ok = —5° 0. This originates from 
(H^^) ~ — 1 MeV. On the other hand, the various phenomenological measurements give 
6k ~ 40° which implies a value of {H~^ ) ~ 40 MeV. Comparing these numbers to {H^^ ) ~ 
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47 MeV extracted from Ref. P] one can see that by extracting a value for {H~^ ) from 9k and 
comparing it to the value for {Hf^ ) one can obtain information about the relative strengths 
of the Coulomb and confining pieces of Hq°p^ . Given the sensitivity of the mixing angle to 
the delicate cancellation between terms, L — S mixing can therefore be a useful means of 
probing the confinement potential.^ 

We next consider the heavy quark limit, where niq oo. In this limit the mass formulae 
simplify to: 

M{^P,) = Mo + {HI J 

M(iPi)/ V V2{HIJ Mo ) VP^ 

Mi'Po)=Mo-2{HlJ (42) 

The two mixed Ki mass eigenstates of J"^ = 1^ appropriate to the heavy quark limit are 
described by the total angular momenta j of the light quark with j = 1/2 and j = 3/2 which 
are degenerate with the = O"*" and = 2+ states respectively. In what follows we will 
take {Hj^) positive but similar results are obtained for {Hj^) negative. For {Hj^) > the 
mixing angle is given by sin6'x = —^2/3 and cos^i^ = y^l/3 {Ok = —54.7") with Mki^^ 
degenerate with the ^Pq (j = 1/2) state and Mk^^^^ degenerate with the ^P2 (j = 3/2) state. 

For the Ki decay constants, in the limit that mg becomes infinitely heavy, the Jk^ become 
proportional to the inverse of the light quark mass and are given by 

fK,.,, = (43) 

So in the heavy quark limit only the j = 1/2 state couples to the weak current. By comparing 
this result to the measured decays one might learn how well the heavy quark limit describes 



^Nevertheless, one cannot rule out the possibility that another mechanism is responsible for ^Pi 
^Pi mixing such as mixing via common decay channels |29|. 
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the strange axial mesons. Using the value of 9k that gives the j = 1/2 and j = 3/2 
eigenstates (expected in the heavy quark limit) and using a finite mass strange quark (still 
taking {H^^ ) > 0) the decay constants are given by: 

Uu = 1/2) = +^(1 + • 



rriu Snis 

/.(.=3/2) = -/fA. (44) 

The value does not change very much for the j = 1/2 state, rriu/Sms ~ 0.2, but the j = 3/2 
state decay constant is no longer zero but is now similar in magnitude to that of the j = 1/2 
state. 

More importantly, the 6^ we used in the above discussion was based on the = 
mass matrix obtained for the heavy quark limit which assumes that the contact and tensor 
contributions are negligible. However, values for these terms extracted from predictions of 
the relativized quark model are: (Hcont) = 33 MeV, {Hten) = 56 MeV, and {H^^^} = 
47 MeV. Clearly the assumption that the contact and tensor pieces are negligible is not 
supported by this model so that the heavy quark limit is questionable for the s quark. 

We conclude that while the heavy quark limit is an interesting means of making qualita- 
tive observations the actual situation for the strange axial mesons is far more complicated. 



VI. CONCLUSIONS 

In this paper we studied the properties of the strange axial mesons in the quark model. 
We extracted the Ki{^Pi)—Ki{^Pi) mixing using the mass predictions, by comparing a quark 
model calculation of the Ki decay constants to the decays r —>■ VtKi and by comparing 
strong decay widths calculated using the pseudoscalar emission model and the fiux-tube 
breaking model to experimental results. In all cases we obtained a mixing angle consistent 
with 9k — 45". There are two important conclusions we can draw from this result. First, 
the relativized quark model predicts a much smaller mixing angle of ~ —5°. Either the 
quark model result is way off, which is possible given the delicate cancellation taking place 
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between the contributions to the spin-orbit term, or a different mechanism is responsible for 
the ^Pi Pi mixing ^9[. The second observation we make on the basis of the quark model 
results is that the heavy quark limit does not appear to be applicable to the strange axial 
mesons. We come to this conclusion because the tensor interaction is still comparable in 
size to the spin-orbit interactions and additionally, the mixing angle is not compatible with 
that expected in the heavy quark limit. 
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FIG. 1. The T ^ Kii' decay widths as a function of Ok for the non-relativistic results. 
R = BR{t vKi{ITIQ))/BR{t i^i^i(1400)). In ah figures the dashed curve is for the wave- 
functions from Ref. [|l8| and the dot-dashed curve is for the wavefunctions from Ref. |17|. In the 
figure for R both curves he on each other. The solid and dotted lines are for the experimental val- 
ues and their 1 — a errors from the TPC/Two-Gamma measurement [Q]. In addition, in the figure 
for BR{t i^Ki), the solid line bounded by the dot-dot-dashed lines are for the CLEO/ALEPH 
result and their 1 — a error [Q. 
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FIG. 3. Predictions of the pseudoscalar emission model as a function of 9k for the Ki partial 
widths (to K*ii and pK), and ratio of D to S amplitudes (to K*!:). The solid curves are for the 
Ki{\27Qi) and the dashed curves are for the iri(1400). The horizontal lines are the 1 — a error 
bounds for the experimental measurements with the same line labelling (solid, dashed) as the 
predictions. (The experimental lower bound for iiri(1400) pK lies on the axis.) 
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FIG. 4. Predictions of the ^Pq model as a function of Ok for the Ki partial widths (to K*'it 
and pK), and ratio of D to S amplitudes (to K*!^). The line labelling is the same as in Fig. 3. 
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TABLES 

TABLE L Axial vector decay constants in the nonrelativistic limit using simple harmonic 
oscillator wavefunctions. The fxi for \K~^) = — \us) are given in units of GeV^. 

Parameter set Mk^ = 1.273 GeV M^i = 1.402 GeV 

Set 1" fK, = -0.049 Ik, = -0.050 

fK. = +0.269 fK. = +0.283 

Set 2^ fK, = -0.070 Ik, = -0.073 

fK. = +0.396 fK^ = +0.415 



From ReL with /? = 0.257 GeV, m„ = = 0.33 GeV and m, = 0.55 GeV. 
^ From ReL |T^ with [3 = 0.3 and = m^^ and as above. 



TABLE IL Axial vector decay constants using the relativized mock-meson matrix elements. 
Simple harmonic oscillator wavefunctions are used with the parameters given below. The fxi for 
\K~^) = — I us) are given in units of GeV^. 

Parameter set Mk^ = 1.273 GeV M^i = 1.402 GeV 

Set 1"^ fK, = -0.024 Jk, = -0.025 

fK. = +0.220 Ik. = +0.231 

Set 2* fK, = -0.040 Ik, = -0.042 

fK^ = +0.486 fK^ = +0.510 



" From ReL ^ with /? = 0.3, m„ = = 0.33 GeV and m, = 0.55 GeV. 
^ We used effective oscillator parameters from ReL . They were obtained by fitting simple 
harmonic oscillator wavefunctions to the rms radii of the wavefunctions of ReL to obtain 
Pk^ = 0.45 GeV. m„ = = 0.22 GeV and = 0.419 GeV. 
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TABLE III. Strong decay amplitudes for the strange axial mesons using the pseudoscalar emis- 
sion model, the ^Pq decay model and the flux-tube breaking model. Note that the amplitudes 
include phase space and are all given in units of MeVi/2_ 



Amplitude 


JT bcULLUoCcildji 


3Pn 

-HD 


Flux-tube breaking 






emission 


Qof 1 0' 
Dfih 1 


Set r 


Set 2'' 


7 








10.4 


12.8 




pK) 


^ 09 
O.UZ 


8 SI 
o.ol 


8.50 


11.0 




pK) 


U.U ( 4 


U.UOD 


0.057 


0.055 




K*Tr) 


lo.o 


io.o 


14.8 


20.5 




K*7r) 


9 9"^ 


9 IP, 
Z.oO 


2.43 


2.29 


g^j^high _ 


>pK) 


15.4 


15.5 


14.8 


20.1 


j^^j^htgh _ 


-pK) 


2.18 


2.28 


2.34 


2.25 


S^j^high _ 


> K*7r) 


17.3 


15.1 


14.2 


20.7 




K*tt) 


4.45 


4.61 


4.74 


4.44 




> Kuj) 


15.1 


15.4 


14.7 


20.0 


jj^j^high _ 


Kuj) 


1.96 


2.04 


2.10 


2.02 


" Simple 


harmonic oscillator wavefunctions 


with (3 = 


0.40 GeV, rriu 


= 0.33 GeV, and 



TJis = 0.55 GeV. 
Wavefunctions from Ref. 
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TABLE IV. Partial decay widths and ratios of D to S amplitudes of the strange axial mesons for 
the pseudoscalar emission model, the ^Pq decay model and the flux-tube breaking model using the 
fitted value of Ok- The widths are given in MeV. Au/As refers to the ratio of D to S amplitudes. 



The errors on 


9k are 


1-a. 












Decay 






Experiment 


Pseudoscalar 




Flux-tube 


breaking" 








(RPP) 


emission 


Set 1^ 


Set l'' 


Set 2° 


9k 








48" ± 5° 


45° ± 4° 


44° ± 4° 


51° ± 3° 


r(Ki(1270) - 


>pK) 




38± 10 


63 


75 


70 


121 


r(Ki(1270) - 


> K*-k) 




14.4 ± 5.5 


16 


12 


11 


35 


lAo/AsiKiimO) ^ 


K*7r)\ 


1.0 ±0.7 


0.64 


0.89 


1.02 


0.40 


r(Ki{im) - 


>pK) 




5.2 ±5.2 


7.9 


12 


13 


6.7 


r{Ki{uoo) - 


> K*-k) 




164 ± 16 


286 


221 


197 


400 


r(is:i(i4oo) - 


> Kuj) 




1.7± 1.7 


2.3 


3.6 


3.9 


1.9 


\Ad/As{Ki{UOO) ^ 


K*-k)\ 


0.04 ±0.01 


0.058 


0.052 


0.051 


0.062 



Note that because the flux-tube breaking calculation involves a numerical integral with 
a 1% error, the two values of S (or D) calculated from the i^i(1270) and -ft'i(1400) decay 
results may not agree exactly. In Table III an average value of the two results is given. 



Because the S and D values in Table III are not exact, using them with Eqn. |2g will not 
exactly reproduce the flux-tube breaking results show in this table, which are calculated 
directly from the numerical work. 

^ Simple harmonic oscillator wavefunctions with /3 = 0.40 GeV, = 0.33 GeV, and = 
0.55 GeV. 

^ Wavefunctions from Ref. fBl. 
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